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Here we report results of a density-matrix-renormalization-group (DMRG) calculation of the 
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well as the rung and leg oxygen atoms are included along with near-neighbor and oxygen-hopping 
matrix elements. This system allows us to study the effects of hole and electron doping on a system 
which is a charge transfer insulator at a filling of one hole per Cu and exhibits power law, d-wave-like 
pairing correlations when doped. In particular, we focus on the differences between doping with 
holes or electrons. 
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I. INTRODUCTION 

Two-leg CuO ladder materials are known to exhibit 
some of the basic physical properties of the high T c 
cuprates. Undoped two-leg ladder materiala-iare found 
to be spin-gapped, charge transfer insulatorsBcl and su- 
perconductivity has observed in hole-doped ladders 
under high pressure Buu Moreover, ladder models have 
proved amenable to numerical studies, allowing one to 
explore the relationship between the parameters in the 
model and physical properties such as the charge and spin 
gaps of the insulating state and the pairing correlations in 
the doped state&Ll Here we continue this type of theoret- 
ical study by examining the properties of a CuO ladder- 
using a density-matrix-renormalization-group (DMRG) 
analysis Ja We will focus on the differences between hole 
and electron doping. 

The geometry of the two- leg ladder which we will study 
is illustrated in Fig. The Cu sites are characterized by 
a d x 2_ y 2 orbital, the rung O sites by a p y orbital, and 
the leg O sites by a p x orbital. The hole Hamiltonian for 
this model is given by 

H = Apd^ptaPiv - t pd \ d ivPi<r +P]a d i°) 

ia (ij) a 

-tpp ^ (pLP]° +P]aPi<r) 

+u d 4^4i d *i + Up^pl-fPitP^Pu (i) 

i i 

Here t p d is the one-hole hopping matrix element between 
nearest-neighbor Cu and O sites and the first sum in 
Eq. (Q) is over all nearest-neighbor sites. We will assume 
for simplicity that t p d has the same value between all the 
Cu and O sites. The second term in Eq. (Q) sums over all 
nearest-neighbor pairs of O sites and t pp is the hopping 
matrix elements between these sites. We have chosen the 




FIG. 1: Schematic lattice structure of a two-leg CuO lad- 
der. Here the solid circles represent Cu(d x 2 _ y i) orbitals and 
the open circles represent 0(p y ) orbitals on the rungs and 
0(p x ) orbitals on both legs. The orbital phases are chosen so 
that there is a hopping matrix element t v d between nearest- 
neighbor Cu and O sites represented by the solid lines and a 
hopping matrix element t pp between nearest-neighbor O sites 
represented by the dashed lines. 



phases of the orbitals such that the signs of the hopping 
matrix elements are constant and with the minus sign 
convention of Eq. (Q), t p d > and t pp > 0. The energy 
difference between the O and Cu sites is A pc i = e p — e d > 
and Ud, and U p are the on-site Coulomb energies for Cu 
and O, respectively. In this model an undoped ladder 
corresponds to a density of one hole per Cu site. We will 
work in units where t p d = 1 and take as typical values 
Ud = 8 and U p = 3, throughout. 

Two of the authors previously studied a similar modelH 
in which the outer O sites were missing and only the 
near-neighbor Cu-0 hopping t p d was present. While this 
simplified form allowed one to compare a charge-gap- 
insulator model with the traditional one-band Hubbard 
and t-J models, the absence of the outer ftpfd the t pp 
hoppings altered the electronic structure. t!JO In addi- 
tion, the electron-doped pairing response of the model 
appeared to be quite different from the hole-doped be- 
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havior. Here, for the more realistic structure shown in 
Fig. |], we return to the study of the two-leg CuO lad- 
der in which the charge-transfer nature of the insulating 
state is treated along with the effects of the oxygen t pp 
hopping matrix elements. We also will study longer lad- 
ders with up to 32 x 2 Cu atoms corresponding to 226 
total O and Cu sites. 

We calculate the static properties-.of the model (|l|) nu- 
merically using the DMRG method!] With this approach 
we obtain accurate ground state energies and expectation 
values (such as correlation functions) for a fixed number 
of holes of each spin. In the DMRG calculations, open- 
end boundary conditions are used so that when we dis- 
cuss local quantities they will actually be averaged over 
the ladder and correlation functions will be calculated 
using distances taken about the midpoint. We have used 
up to to = 2000 density matrix eigenstates to build the 
DMRG basis. Using an extrapolation of DMRG ground 
state energies to vanishing discarded weight ,E3 we have 
obtained ground state energies and gaps which are accu- 
rate to parts in 10~ 3 t p d- Although the largest source of 
errors in our calculations are finite size effects, we have 
been able to extrapolate a number of gap measurements 
to an infinite-size ladder. 

In Sec. [n] we discuss our results for a filling of one hole 
per Cu, which corresponds to the undoped ladder. Our 
main interest is in the effect of the site energy difference 
A p d = e p — e<i and the oxygen hopping parameter t pp . 
We will examine the charge and spin magnitudes on the 
various sites and the dependence of the charge and spin 
and t pp . When Ud is large compared with 



> 2t pp , we find as expected that the ladder 



gaps on A p 
A pd and A p 

is a charge transfer insulator with a spin gap. By com- 
paring the low-lying spin states of the CuO model with a 
two-leg Heisenberg model, we extract effective rung and 
leg exchange interactions. In Sec. [II, we examine the 



charge and spin magnitude distributions for the doped 
case and find that, for typical parameters characteristic 
of the cuprates, the doped holes tend to be spread out on 
the O lattice while doped electrons tend to be localized 
on the Cu sites. We also study the effect of doping on the 
effective exchange interactions and on the local spin-spin 
correlations. We then turn to a study of the pairing cor- 
relations for both the hole-doped and the electron-doped 
ladder. We find a power law behavior with negative "d- 
wave-like", Cu rung-leg pair correlations for both hole 
and electron doping. However, we find that the internal 
structure of a pair depends upon whether the system is 
electron- or hole-doped. In Sec. IV, we conclude with 
some further comments relating these results to the gen- 
eral problem posed by the cuprate materials. 



II. THE UNDOPED LADDER 

We first investigate the properties of an undoped lad- 
der. For local properties such as the charge and spin 
magnitude distributions on the sites and the spin corre- 
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FIG. 2: Hole density (ncu) and (no) on the various sites 
of an undoped ladder with 8 x 2 Cu atoms: Cu sites (solid 
circles), rung 0(p y ) sites (up-triangles) , leg 0(p x ) sites (down- 
triangles) , and outer O sites (squares) versus (a) t pp for A p d = 
2.7 and (b) A p d for t pp = 0.5. 



lations between sites, a ladder containing 8x2 Cu sites 
is sufficiently long. In Fig. |^ we show the average ground 
state hole occupation of Cu sites (ncu) — (d\^ d^ + d\^dii) 
and of the various O sites (no) = (pjtPit + PiiPil) as a 
function of t pp and A 
tion is shown for A p d = 2.7 versus t pp 
the hole occupation on the Cu sites decreases and that 
on the O sites increases. In Fig. ^(b), the hole occu- 
pation is shown for t pp — 0.5 versus A p d- Here as A p d 
increases, the hole occupation on the Cu sites increases 
as one would expect. The hole occupation on the outer 
O sites is typically about half that on the other O sites. 
For tpp — 0.5 and A p d — 3, the holes are approximately 
70% on the Cu sites and 30% on the O sites. 

The average value of the square of the spin moment 



pp aim L-L P d- In Fig. |2|(a) the hole occupa- 

As t pp increases 



(2) 



on various sites is shown in Fig. |^ for the same parame- 



ters. Here = d] a or p\ a depending upon the site, <x is 
the usual Pauli spin matrix, and the indices a and /3 are 
summed over. For A p d = 3 and t pp = 0.5, the average 
spin moment on the Cu site is 0.5. Dividing this by the 
average hole occupation for the Cu gives a squared spin 
moment of order 0.7 compared with 3/4 for a Heisenberg 
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FIG. 3: Local spin moment (S^) on the various types of 
sites of an undoped ladder with 8x2 Cu atoms: Cu sites 
(solid circles), rung 0(p y ) sites (up-triangles) , leg 0^) sites 
(down-triangles), and outer O sites (squares) versus (a) t pp 
for A p d = 2.7 and (b) A p d for t pp = 0.5. 



spin s = 1/2. With Ud — 8, the spin moment is well 
formed on the Cu site and (S^ u ) simply tracks (nc u ); 
the probability of having a hole on a Cu site. 

The charge gap of an undoped ladder with L x 2 Cu 
sites is given by 

A C (L) = E (+2, L) + E (-2, L) - 2£ (0, L) , (3) 

where E (N, L) denotes the ground state energy of a lad- 
der of length L with N + 2L holes. Extrapolating the 
charge gap from numerical data for up to L = 16, we 
find the results for A C (L — » oo) plotted in Fig. ||. Here 
Fig. ||(a) shows A C (L — > oo) versus t pp for various values 
of A p d and Fig. ^(b) shows A C (L — » oo) versus A p d for 
various values of t pp . From these results, we see that in 
order to hav£_a_charge gap, we need A p d > 2t pp . Ex- 
perimentallyBEJ, A c is of order 1.4 to 2 eV. If we take 
tpp ~ tpd/2 thenj-Apd w 3 gives a realistic charge gap 
with t pd w 1.6eV.E3 

The lowest triplet spin excitation is gapped and local- 
ized on the Cu sites of the ladder (i.e., the spin density 
remains zero on the O sites in the triplet state). The spin 
gap is given by 

» , fN N \ f N N \ 

A s (N,L) = E ( K - + l,--l,Ly Eo ^,-,L^ 

(4) 




FIG. 4: Charge gap of undoped ladders extrapolated to an 
infinite ladder length L — > oo. (a) As a function of t pp for 
A p d = (squares), A p d = 1 (down triangles), A p d = 2.5 (up 
triangles) , and A p d = 4 (circles) . (b) As a function of A p d for 
t PP = (circles), t pp — 0.5 (triangles), and t pp = 1 (squares). 



where Eo(Nf , N±, L) is the ground state energy of an L x 2 
ladder with N a + L holes of spin a =|, |. We first discuss 
the spin gap of the undoped ladder (JV = 0). In this case 
we extrapolate the spin gap to an infinite ladder length 
using data for up to L — 16. The dependence of the spin 
gap A S (L — > oo) on £ pp and A p d is shown in Figs. |](a) 
and ||(b), respectively. As for the charge gap, there is 
a spin gap only if A p d ^ 2t pp . For large A p d the spin 
gap decreases due to a decrease in the effective exchange 
coupling as observed in our previous work.13 Thus the spin 
gap goes through a maximum at a finite value of A p d > 
2t pd . Experimentally,BEl A s - 0. 03-0. 05eV so that with 
t pd ~ 1.6eV, we need A s /t p d — 0.03. This confirms that 
t pp ~ 0.5 and A p d ~ 3 represent appropriate choiccS|, 
consistent with electronic band structure calculations £3 

For a two-leg Heisenberg ladder, the spin gap is equal 
to approximately half of the exchange interaction. Thus 
from Fig. @(a) we see that provided A p d ^ 2t pp , the 
effective exchange interaction between Cu spins increases 
with t pp . For example for A p d = 2.5, the spin gap and 
hence the effective Cu-Cu exchange increase by 40 % as 
tpp goes from to 0.5. This is due to the additional 
exchange path through the O ions. This crucial role of 
the 0-0 hopping t pp in giving a large effective Cu-Cu. 
exchange has been emphasized by Eskes and Jefferson.Ej 
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FIG. 5: Spin gap of undoped ladders extrapolated to an 
infinite ladder length L —* oo. (a) As a function of t pp for 
A p d = 1 (circles), A p d = 2.5 (up triangles), A p d = 4 (down 
triangles), and A pt j = 7 (squares), (b) As a function of A ye j 
for t pp = (squares), f pp = 0.5 (up triangles), and t pp = 1 
(circles). 



The nearest-neighbor Cu spin-spin correlations (Si-Sj) 
are antiferromagnetic in the undoped ladder. For a fixed 
value of tpp, —(Si ■ Sj) increases with A p d and tends to 
the value obtained for a two-leg Heisenberg ladder when 
A p d ^> tp P . This simply reflects the fact that increasing 
Apd causes the holes to become more localized on the Cu 
sites, as shown by the increase of the Cu spin moments 
(Sf) in Fig. ||(b). For the same reason, for a fixed value of 
Apd, — (Si-Sj) decreases with increasing t pp and vanishes 
for tpp 3> Apd- This decrease is due to the reduction of 
the local magnetic moments (Sf) on the Cu sites as the 
holes spread further onto the O sites, as seen in Fig. ||(a). 

The low energy properties of the spin degrees of free- 
dom of the undoped d-p Hubbard model (|I|) can be 
mapped onto an anisotropic Heisenberg ladder. 

H = J2Jii Si-Sj. (5) 

Here (ij) represent near-neighbor bonds along the legs 
and rungs of a 2-leg ladder with Jy = J± for the rungs 
and Ji| for the legs. To estimate the parameters J± 
and Jy , we consider the low-energy states of the d-p 
model for a CU2LO5L cluster (L x 2 Cu ladder). Then 
by making a one-to-one correspondence of the states of 
this system to the L x 2 Heisenberg ladder, we derive 



the Jl and Jy parameter of the anisotropic Heisenberg 
Hamiltonian. These processes are performed for ladders 
of length L = 8, 12, and 16 with the open boundary 
condition. We confirm that the low-lying spin states 
of the d-p and Heisenberg clusters have the same quan- 
tum numbers and characters for all the three values of 
L. To make the mapping systematically with various 
system sizes, we extract the two lowest triplet excita- 
tion energies, E tl = E (S to t = 1) - E (S tot = 0) and 
E n = Ei (Stot = 1) - MS = 0). Here, E n (S = S tot ) 
is the energy of the n th excited state with total spin S. 
Then we determine J± and J|| by fitting the excitation 
energies E t \ and £"42 with results obtained from d-p clus- 
ters with the same L. In this way, with A p d = 3 and 
tpp = 0.5, we find that J± = 0.06i p( z(~ 96meV) and 
J11 = O.QStpdi^ 128meV), both of which are in a reason- 
able ra nae as. compared to experimental and theoretical 
resultsJilHy 

It is also interesting to examine the effect of t pp and 
Apd on the effective exchange coupling anisotropy. From 
the nearest-neighbor Cu spin-spin correlation (Si ■ Sj) we 
define 

(Si ■ S,) 

„ \ ' J ' rung /„x 

R ~ (s.-s^^ (6) 

where we use an average of (Si ■ Sj) over the ladder. 
For a Heisenberg ladder (|B|) with an exchange coupling 
Ji on the rungs and Ju on the legs, R versus is 
shown in Fi g . J q(a) . Various experimental and theoreti- 
cal estimatea3oE3 suggest that values of J± / J|| ranging 
from 0.5 to 1.0 are appropriate for two-leg CuO ladders. 
Figure ||(a) shows that this corresponds to values of R 
varying from 0.6 to 1.2. The ratio R obtained for our 
CuO ladder are shown versus t pp and A p d in Figs. ||(b) 
and ||(c), respectively. We see that the exchange cou- 
pling anisotropy is again consistent with parameter val- 
ues Apd ~ 3 and t pp ps 0.5. In the absence of the 
outer O sites, the anisotropy ratio R becomes signifi- 
cantly larger than the values shown in Figs, ^(b) and 
^(c) for < tpp < 1. In the charge transfer insulat- 
ing phase, one finds an effective exchange coupling ratio 
J±/J\\ > 1 for realistic parameters, in disagreement with 
most experimental and theoretical estimates. Thus, it 
appears that it is important to include the outer O sites 
in order to obtain a qualitatively correct description of 
two- leg CuO ladders. 

The dependence of the charge and spin magnitude dis- 
tributions, the charge and spin gaps, and the spin-spin 
correlations on the parameters t pp and A p d can be under- 
stood qualitatively. For t pp > A p d/2, the holes are delo- 
calized on all Cu and O sites and the undoped CuO ladder 
is in a metallic phase. For A p d ^ 2t pp , the holes tend to 
be localized on the Cu sites and the undoped system is 
a charge-transfer insulator. Low-energy spin excitations 
involve the holes on the Cu sites only. Thus, the spin 
degrees of freedom localized on the Cu sites lead to an 
effective two-leg spin ladder. Here, as noted, we are work- 
ing in units where t p d = 1, Ud = 8, and U p — 3. We have 
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FIG. 6: Ratio R between nearest-neighbor spin correlations 
on rungs and legs, (a) For a two-leg Heisenberg ladder with 
100 rungs as a function of J± / Jj| . For an undoped CuO ladder 
with 8x2 Cu sites (b) versus t pp and (c) versus A p( j. 

seen that the magnitudes of the charge and spin gaps as 
well as the site charge densities, spin moments, and ef- 
fective exchange couplings are consistent with A p d = 3 
and t pp = 0.5. We will now examine what happens when 
such a CuO ladder is doped with holes and electrons. 



III. THE DOPED LADDER 

Turning to the doped situation, we investigate the 
properties of the model ([!]) with A p d = 3, Ud — 8, U p = 3, 
and tpp = 0.5 for various hole concentrations (per Cu 
atom) x = 1 + N/(2L), where N is the number of doped 
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FIG. 7: (a) The hole density n and (b) the square of the spin 
moment (Sf) on the various types of sites of a ladder with 
16 x 2 Cu atoms as a function of filling x : Cu sites (solid 
circles), rung 0(p y ) sites (up triangles), leg 0(p x ) sites (down 
triangles), and outer O sites (squares). 



holes (N > 0,x > 1) or doped electrons (N < 0, x < 1) 
in a ladder with L x 2 Cu atoms. Results for the site 
occupations versus the doping are shown in Fig. f7](a). 
Here, we see that the slope of the hole occupation on the 
Cu site versus x changes at x = 1. For electron doping, 
the electrons are more likely to go on the Cu site than 
are the holes when the system is hole doped. This is, of 
course, what one would expect for a charge-transfer insu- 
lator. This is also consistcnt-with Monte Carlo results for 
the 3-band Hubbard modelOta as_well as our previous 
DMRG results for the CuO ladder The variation of the 
local squared spin moments on the different sites versus 
doping are shown in Fig. 0(b). As discussed previously, 
in the charge transfer regime, where Ud is large compared 
to Apd and A p d is larger than several times t pp , the spin 
moments are mainly on the Cu sites. When the system is 
electron doped, the electrons go onto the Cu sites elim- 
inating their moments. Thus, the average moment on 
the Cu sites depends upon x in the same way as the Cu 
site occupation (rtcu), and (Sq u )/{71cu) remains essen- 
tially constant, of order 0.7. The change in the average 
moment shown in Fig. 0(b) simply reflects the fact that 
for the electron-doped case, spin moments at individual 
Cu sites are removed when electrons are added to those 
sites. Likewise, for the hole-doped situation the average 
spin moment on the Cu sites increases only slightly as x 
increases and a few more Cu sites have holes. 
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FIG. 8: The average nearest-neighbor Cu spin-spin corre- 
lation function {Si ■ Sj) on the rungs (solid circles) and legs 
(open circles) versus x calculated on ladders with 16 x 2 Cu 
sites. 



In Fig. |^ we have plotted the average value of the 
nearest- neighbor Cu spin-spin correlations (Sj • Sj) on 
the rungs and legs versus doping. As electrons are doped 
onto the ladder and spin moments are removed from the 
Cu sites, the magnitude of the (Si ■ Sj) correlations de- 
crease on both the rungs and the legs. Similar results 
are observed for the nearest-neighbor spin-spin correla- 
tions in the t — J model on a two-leg ladder as the hole 
concentration increases. 

However, as holes are added to the CuO Hubbard 
model ([!]) and x increases, the rung spin-spin correla- 
tions decrease more rapidly than those along the legs. 
As we have seen, the,-added holes tend to go onto the O 
sites which frustrates^ or screens the exchange coupling 
between the Cu spins. A rough estimate of this can be 
seen by considering a CU4O12 cluster with one hole or 
electron added. Calculating the low-lying states of the 
cluster and comparing the excitation energies with a 4- 
site t-J e g system we find that J e s/t p d is 50% smaller for 
the one hole (x = 1.25) case relative to its value for the 
one electron (x — 0.75) case. Thus, in the hole-doped sys- 
tem, the effective exchange interaction between the spin 
moments on the Cu sites is weakened. Because this effec- 
tive exchange is initially larger along the legs than across 
the rungs, the rung correlations decrease more rapidly 
than those along the legs. We note that the anisotropy 
in the reduction of the spin-spin correlations as a func- 
tion of doping is much stronger than in a t — J model on 
a two- leg ladder, even if one uses a strongly anisotropic 
exchange coupling J± ^ J». 

The average values of the nearest-neighbor Cu-0 spin- 
spin correlations increase rapidly with the number of 
holes in the ladder. For hole doping the appearance 
of strong nearest-neighbor Cu-0 spin correlation&.corre- 
sponds to the formation of a Zhang-Rice singletEll in an 
isotropic two-dimensional Cu02 lattice. In the two-leg 
CuO ladder, however, the strength of the Cu-0 spin-spin 
correlation depends on the type of the O site considered. 
Rung Cu-0 spin-spin correlations increase less rapidly 
than the other ones with increasing hole concentration 




FIG. 9: The Cu rung-rung pair held correlation function 
D{£) = (A i+ iAl) versus I for a 32 x 2 Cu ladder. Here, 
we have normalized D{£) to its value at £ = 4. Circles, 
squares, and triangles correspond to an undoped ladder, a 
ladder doped with two holes, and with two electrons, respec- 
tively. Inset: Same data (normalized to their £ = 1 value) on 
a log-log scale. The lines have slope —1/2, —1 and —2. 

x. This confirms the higher frustration of antiferromag- 
netic correlations in the rung direction than in the leg 
direction when an additional hole is placed between the 
two holes localized on nearest-neighbor Cu sites. Cor- 
respondingly, for hole doping the spin-spin correlation 
between a Cu site and its outer O site increases faster 
than for the other nearest-neighbor sites because of the 
absence of frustration for this type of O sites. 

Turning next to the pairing correlations, we have cal- 
culated the Cu rung-rung pair field correlation function 

D(£) = (A t+e A\) (7) 

with 



Here, d iXs creates an electron of spin s on the i rung 
and the A = 1 or 2 leg of the ladder. The pair held has ll d- 
wave-like" structure in the sense that the rung-leg Cu-Cu 
pair held correlation function, has a negative sign while 
the rung-rung or leg-leg pairing correlations are positive. 
We have also examined other pair held correlation func- 
tions corresponding to pairs on diagonal Cu sites or pairs 
on O sites and found that they are either qualitatively 
similar to the Cu rung-rung pair field correlations or de- 
crease much faster. Some results for the rung-rung pair 
field correlation function D{£) versus I are shown in Fig. ^ 
for a ladder containing 32 x 2 Cu sites. We have normal- 
ized D{£) with respect to its value at £ = 4. The pair 
field correlations in the undoped system decay rapidly 
(exponentially) with distance while, as shown in the in- 
set, the pair correlations when two electrons are added or 
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removed exhibit an approximate power law decay. (Qual- 
itatively similar pair field correlation functions have been 
calculated in a t— J two-leg ladder with two doped holes.) 
While the decay of the normalized rung-rung pair field 
correlations for the 2-electron and 2-hole doped ladders 
are similar in the CuO Hubbard model (0), the size of 
the rung-rung Cu pairing correlations are a factor of or- 
der 4 larger for the electron-doped case relative to the 
hole-doped one. This reflects the fact that for electron 
doping, the added carriers go primarily onto the Cu sites 
while for hole doping they go primarily onto the O sites. 
For higher doping (0.75 < x < 1.25) the Cu rung-rung 
pair field correlations are very similar to those shown in 
Fig. H for doped holes or electrons. 

To determine whether doped holes or electrons form 
bound pairs we have first calculated the binding energy of 
two doped holes and electrons. The pair binding energy 
is defined as 

A pb (L) = 2E (±1, L) - E {±2, L) - E (0, L), (9) 

where Eo(N,L) is the ground state energy of a ladder 
with L x 2 Cu site and N + 2L holes. We could not 
extrapolate A p b(L) to the L — > oo limit because finite- 
size effects are too irregular. However, we have found 
that A p b(L — > oo) is certainly positive both for hole 
and electron doping. Therefore, we expect two doped 
holes or electrons to build a bound pair in the CuO Hub- 
bard ladder for_the parameters we use here. As in our 
previous work,fj the binding energy of two doped holes 
is of the order of the spin gap of the undoped ladder. 
For instance, for a 16 x 2 Cu site ladder we have ob- 
tained A p b rj O.Q2t p d, to be compared with the spin gap 
A S (L — > oo) = O.QStpd- The binding energy of two elec- 
trons, Apt w 0.04t p( j, is also of the order of the spin gap 
in the undoped ladder but twice as large as the hole pair 
binding energyr-Jn our previous investigation of the two- 
leg CuO ladderawe found a vanishing pair binding energy 
for electron doping. To understand this change we have 
calculated the pair binding energy as a function of the 
various model parameters in a ladder with 8x2 Cu sites. 
It appears that the difference between the present result 
and the previous oneEl is due to the nearest-neighbor O 
hopping tpp. In fact, t pp > is necessary to obtain a 
finite pair binding energy in the electron doped ladder. 
When the CuO model is reduced to a simple, one-band 
Hubbard model, t pp > leads to an effective next-near- 
neighbor hopping which favors pairing when the system 
is electron doped. c3 In addition, as previously discussed, 
tpp > leads to an enhanced Cu-Cu effective exchange. 

In order to have a more detailed picture of the differ- 
ence between the electron- and hole-doped pairs, we have 
investigated their internal structure. In a t — J model or 
a spin-fcrmion model, it is easy to measure correlations 
between doped charges because there are no other itin- 
erant charges (i.e., there are no charge fluctuations in 
the undoped system). In the CuO Hubbard model ([!]), 
however, bare electron-electron and hole-hole correlation 
functions give little information about the correlations 



t — I I I | I I I I | I I I I 




5 10 15 



FIG. 10: The projected correlation functions C e ,ji(0 for (a) 
the Cu sites of a two-electron doped ladder and (b) the O4 
orbitals of a two-hole doped ladder. The distance I is mea- 
sured along the ladder in units of the interval between two 
Cu rungs. Open and solid circles correspond to correlations 
between doped charges on the same leg and on different legs, 
respectively. 

between two doped charges because they are dominated 
by the quantum fluctuations of the other charges. There- 
fore, we have to project out the ground state of (|l|) onto 
a subspace with no charge fluctuation. For that purpose 
we consider a perturbation expansion in t pc i and t pp of 
the Hamiltonian (Q) with Ud > A p d > and x close to 
1 . The ground state of the unperturbated Hamiltonian 
(tpd = tpp = 0) is degenerate. For electron doping the 
ground state has at most one hole on each Cu site and 
no hole on any O site. The projection operator onto the 
ground-state subspace is 

p e = jj(i - 4 T *T<*i) II ;/v/<>,./<;.; ■ (10) 

i i 

For hole doping there is exactly one hole on each Cu 
site and no doubly occupied O site in the ground state of 
the unperturbated Hamiltonian. The projection operator 
onto the ground-state subspace is 

p h = n(s? T + st^ n(i - p} T p 3 -TP im) > (11) 

i 3 
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where Sf a = d\ a d ia di_ a d\ _ a . Using perturbation theory 
one could then derive effective Hamiltonians (a general- 
ized t — J model and a generalized spin-fermion model, 
respectively) in the subspace defined by P e and Pf,, , which 
approximately describe the low-energy properties of the 
CuO model (|f|) in the regime t p d,t pp <C A p d < Ud for 
electron doping and hole doping, respectively. There- 
fore, for the realistic parameters A^d — 3tpd = 6i pp , we 



ft 



(a) 



calculate the ground state \ip) of (|lj) and then use the 
projected states P e \ip) and Ph\ip) to determine the corre- 
lations between doped carriers (electrons or holes). 

Figure 10(a) shows the electron-electron correlation 
function 



C?(m) 



\P e E n E n+m P e \j>) 
(iP\P e E n P e \iP) 



(12) 



where E n — dn^d'^dn^d^,, calculated in a ladder with 
32 x 2 Cu atoms and two doped electrons. C™ (m) mea- 
sures the correlation between two added electrons on the 
Cu sites. Note that the most probable arrangement for 
the two electrons is on diagonal nearest-neighbor sites, 
which is similar to what is seen in t-J ladders. Results 
for the hole-doped case are shown in Fig. |l^(b). Here the 
hole-hole correlation function is given by 



C£(m) = 



(^\P h £ CT H na H n+m ^ a P h \ilj) 



(ip\P h E„ H na P h \iP) 



where the operator 



(13) 



(14) 



(nj) 



measures the probability of finding a hole in an O4 orbital 
around the Cu site with index n (the sum runs over the 
four nearest-neighbor O sites of this Cu site) . We average 
the density over four O sites around each Cu site because 
a doped hole is locked in a singlet state with another hole 
on a C-iLsite as in a CUO4 cluster (i.e., in a Zhang-Rice 
singletEJ). This also facilitates the comparison with the 
electron-doped case and the t — J two- leg ladder. C%(m) 
measures the correlation between two added holes in O4 
orbitals around different Cu sites. Figure |lC](b) shows 
that the most probable arrangement for the holes is on 
diagonal next-nearest-neighbor O4 orbitals. 

For larger distance m, both C™(m) and C£(m) de- 
crease approximately as exp(— |m|/£) with correlation 
lengths (in units of the distance between two Cu rungs) 
£ e w 3.3 and £/, w 2.9, respectively. Note that these 
correlation lengths cannot be interpreted as the pair size 
because the exponential decay sets in only for m £ e ,fr,- 
The average distance between doped carriers is given by 



m e< h 



Em \m\C e ,h(m) 



(15) 



which yields m e w 4.5 and « 4.8 for the electron and 
hole pair, respectively. 




£ 11.07 



FIG. 11: A schematic view of the density-density correlations 
showing the structure of a pair in (a) an electron-doped ladder 
and (b) a hole-doped ladder. Here, one carrier (doped electron 
or hole) is located at the site n marked by an open circle 
(around the middle of the ladder). The radius of the solid 
circles is proportional to the probability C"(m) [C£(m)] of 
finding the second doped electron (hole) on that site. Only 
the central part of ladders with 32 x 2 Cu sites is shown. 



A more intuitive picture of the difference in the internal 
pair structure between electron and hole doping is shown 
in Fig. [ll]. Here, we show the probability C e (m) [Ch(m)] 
of finding the second doped electron (hole) at a Cu site 
(in a O4 orbital) when the first doped electron (hole) 
is located on the Cu site (in the O4 orbital) marked by 
an open circle. The electron pair appears to be denser 
than the hole pair which is qualitatively consistent with 
the larger binding energy (||) calculated for the electron- 
doped case. 

As we use open boundary conditions, the correlation 
functions C™ h (m) depend on the position (i.e., the index 
n) of the first doped carrier (electron or hole) . We have 
checked the pair structure for different ladder sizes and 
at different positions on the ladder and have found sig- 
nificant variations close to the ladder ends. In Fig. [ll] 
we show the pair structure in the middle of a ladder with 
32 x 2 Cu sites as finite-size and boundary effects are min- 
imal there. We believe this structure to be representative 
of the pair structure in an infinite ladder. 

As in the undoped ladder, the lowest triplet excitation 
in the doped ladder is gapped and involves only holes 
on the Cu sites. The spingap (||) for the doped system 
versus x is shown in Fig. H2. Here, we have extrapolated 
the spin gap A S (L) to L — > 00 for fixed hole concentra- 
tions x = 0, ±1/8, and ±1/4 using numerical results for 
up to L = 32. For the limit x — > 1 we have extrapolated 
the spin gap A S (L) to L — » 00 for the 2-hole and the 2- 
electron doped ladders with up to L = 32. One sees that 
the value of the spin gap in the i -> 1 limit differs from 
the x — 1 undoped value. That is, the q = (tt, 7t) magnon 
in the undoped system (x — 1) has a gap of order 0.03, 
while in the doped system the spin gap as x — > 1 is of 
order 0.018. This behavior is similar to wbat has been 
found previously in studies of t-t'-J laddersca. 

An examination of the spin structure of the S z = 1 
state containing two doped electrons or holes shows that 
the spin and added carriers are spatially correlated. 
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FIG. 12: The L — > oo extrapolated spin gap A s versus hole 
concentration x. The open squares denote hole doping while 
the triangles denote electron doping. The undoped x = 1 case 
is shown as the open circle. 



Thus, the spin gap is set by the energy of a bound 
magnon-carrier pair. For instance, Figure [l^ shows the 
spin density (on the Cu sites) of the S z — 1 state with 
two doped carriers (electrons or holes), when these doped 
carriers are in the configuration with the highest prob- 
ability Ce^him) (the pair structure of the triplet states 
is similar to that of the singlet states shown in Fig. 11). 
For electron doping, the spin structure shown in FisJL3|is 
very similar to that found in t-t'-J two- leg ladders£3For 
hole doping, however, the spin structure is quite different, 
which is only in part due to the difference in the carrier 
configuration. Even if we use the same carrier configura- 
tion as for the electron-doped case (i.e., two doped carri- 
ers on nearest- neighbor diagonal sites), the spin structure 
of the hole-doped system is different from that found in 
the electron-doped system or t-J model. In particular, 
we note that there is a small spin density on the Cu sites 
in the center of the O4 orbitals where the holes are lo- 
cated, which shows that in the magnon-carrier pair, holes 
on Cu and O sites are no longer completely locked in a 
Zhang- Rice singlet. 



(a) 



If 0.0 



(b) 



f t i 



\f 0.0 



FIG. 13: A schematic view of a magnon in (a) an electron- 
doped ladder and (b) a hole-doped ladder. Here, both doped 
carriers (electrons or holes) are located at the positions (Cu 
sites or O4 orbitals) marked by open circles around the middle 
of a ladder with 16 x 2 and 25 x 2 Cu sites, respectively. The 
length of the arrows is proportional to the spin density on the 
corresponding Cu site. 



IV. CONCLUSION 

Using the density matrix renormalization group, we 
have studied the structure of the charge, spin, and pair- 
ing correlations for a Hubbard-like model of a two-leg 
CuO ladder. This model allowed us to examine the dif- 
ferences between electron and hole doping that occur in 
a charge transfer insulator. For the undoped ladder, 
where there is one hole per Cu, we found that pa r am- , 
eters obtained from LDA downfolding calculationst2ltLil 
lead to reasonable charge and spin gaps. A typical pa- 
rameter set in units of the Cu-0 hopping t p d were a near- 
neighbor oxygen-oxygen hopping t pp = 0.5, an oxygen- 
copper site energy difference A p d = 3, and on-site Cu and 
O Coulomb repulsions Ud — 8 and U p — 3 respectively. 
With these parameters, the charge gap is determined pri- 
marily by A p d rather than Ud and the spin gap is set by 
the effective Cu-Cu exchange. We have found that the 
0-0 hopping t pp plays a significant role in giving a large 
effective Cu-Cu exchange interaction. 

For the undoped ladder, we found that the holes were 
distributed approximately 70% on the Cu sites and 30% 
on the O sites. The large on-site Coulomb interaction on 
the Cu site leads to the spin moment being dominantly 
on the Cu site. Basically, if the hole is on a Cu site, the 
square of the spin moment (0.7) has nearly its full value 
of 3/4. When electrons are added, they go primarily 
onto the Cu sites (of order 80% of the added charge goes 
onto the Cu sites). This is seen in both the decrease in 
the average Cu site hole occupation and the decrease in 
the average spin moment. Note that the decrease in the 
spin moment simply reflects the fact that there are fewer 
Cu ++ sites in the electron-doped system. Alternatively, 
when holes are added (electrons are removed), they go 
primarily onto the O sites (of order 20% of the holes 
go onto the Cu with 80% going onto the surrounding O 
sites). This leaves the local spin moments on the Cu 
sites. 

Thus, the undoped CuO ladder is a spin-gapped in- 
sulator with a charge gap set by the oxygen-copper site 
energy difference. Magnetically the undoped ladder is 
essentially a Heisenberg ladder made up of hole spins 
localized on the Cu sites. When the ladder is electron 
doped, some of the holes (and spins) on the Cu sites 
are removed and pairing correlations develop. Because 
the doped electrons go primarily onto the Cu sites, the 
electron-doped system is closer to the one-band Hubbard 
or t-J models. When the ladder is hole doped, the holes 
go dominantly onto the O sites leading to a low-density 
gas of fermions delocalized over the O-site lattice and in- 
teracting by spin-exchange with the localized spins on the 
Cu sites. The local magnetic moments remain essentially 
unchanged on the Cu sites but the O holes frustrate the 
exchange coupling between the spins on the Cu sites. 

This local, strong coupling picture focuses on the dif- 
ferences between the electron and hole-doped ladders. 
Nevertheless, on low energy scales, the spin gap and pair- 
ing correlations of the electron- and hole-doped systems 
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are quite similar. Both exhibit d-wave like power law 
pairing correlations in which the pair field Cu-Cu rung- 
leg correlations are negative. The spin gap for both dop- 
ings is associated with a bound magnon-pair. However, 
the internal structure of the pairs differ for electron and 
hole doping. 
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